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Abstract 

We  discuss  the  Navier-Stokes  equations  for  an  incompressible  fluid  with  a  viscosity 
that  is  allowed  to  depend  on  the  pressure.  Ellipticity  and  the  complementing  condition  of 
Agmon,  Douglis  and  Nirenberg  [lj  are  discussed.  It  is  found  that  the  pressure  dependence 
of  viscosity  leads  to  the  possibility  of  a  change  of  type.  It  is  shown  that  the  Dirichlet 
initial-boundary  value  problem  is  well  posed  as  long  as  the  equations  do  not  change  type. 
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Significance  and  Explanation 


In  most  mathematical  treatments  of  the  Navier-Stokes  equations,  it  is  assumed  that 
the  viscosity  is  a  constant.  Viscosities  of  real  fluids,  however,  depend  not  only  on  the  tem¬ 
perature,  but  may  also  change  significantly  with  pressure,  in  particular  at  high  pressures. 
In  this  paper,  the  mathematical  consequences  of  such  a  pressure  dependence  are  investi¬ 
gated.  It  is  found  that,  in  contrast  to  the  ordinary  Navier-Stokes  equations,  ellipticity  can 
be  lost,  and  the  equations  are  not  necessarily  well-posed.  The  complementing  condition 
for  traction  boundary  conditions  is  investigated,  and  an  existence  theorem  for  the  initial¬ 
boundary  value  problem  with  prescribed  velocities  on  the  wall  is  proved.  One  of  the  main 
differences  to  ordinary  Navier-Stokes  theory  lies  in  the  elimination  of  the  pressure,  which 
now  leads  to  a  nonlinear  elliptic  partial  differential  equation  instead  of  Laplace’s  equation. 
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SOME  REMARKS  ON  THE  NAVIER-STOKES  EQUATIONS 


WITH  A  PRESSURE-DEPENDENT  VISCOSITY 

Michael  Renardy 

1.  Introduction 

In  most  mathematical  treatments  of  the  Navier-Stokes  equations,  it  is  assumed  that 
the  viscosity  of  the  fluid  is  constant.  Viscosities  of  real  fluids,  however,  depend  not  only 
on  the  temperature,  but  may  also  depend  on  the  pressure.  In  fact,  the  possibility  of 
a  pressure-dependent  viscosity,  even  for  incompressible  fluids,  is  already  mentioned  by 
Stokes  [8].  Stokes  does  not  pursue  the  matter  further,  since  experiments  at  his  time  did 
not  suggest  that  this  pressure  dependence  was  important.  More  recent  experiments  at  very 
high  pressures,  however,  have  shown  considerable  changes  in  viscosity.  In  [3],  for  example, 
various  organic  liquids  were  studied  at  pressures  of  up  to  several  thousand  atmospheres.  It 
was  found  that  for  some  of  these  liquids  the  viscosity  increased  by  two  orders  of  magnitude 
compared  to  its  value  at  atmospheric  pressure,  while  compression  was  only  of  the  order  of 
10%  [2].  Such  data  suggest  that  it  is  reasonable  to  consider  a  pressure  dependence  of  the 
viscosity,  even  in  fluids  that  can,  to  a  good  approximation,  be  regarded  as  incompressible. 
The  results  of  [3}  indicate  that  the  dependence  of  viscosity  on  pressure  becomes  stronger 
as  the  pressure  increases  and  that  it  is  stronger  in  fluids  that  already  have  a  high  viscosity 
at  atmospheric  pressure.  We  note  that  the  viscosities  of  polymer  melts  are  several  orders 
of  magnitude  higher  than  that  of  any  of  the  fluids  studied  in  [3]. 

From  a  theoretical  point  of  view,  we  can  think  of  incompressible  fluids  as  a  limiting 
case  of  compressible  fluids.  For  a  Newtonian  compressible  fluid,  the  stress  tensor  is  given 
by 

T  =  -<£(p)l  +  2>?(p)(D  -  ^  div  u  1)  +  ?(p)  div  u  1.  (1) 

Here  D  denotes  the  symmetric  part  of  the  velocity  gradient,  and  D  -  |  div  u  1  is  the 
deviatoric  (traceless)  part  of  D.  The  total  “pressure”  is  given  by 

P  =  ~  ?(p)  div  u.  (2) 

Noting  that  div  u  =  -  In  p,  we  see  that  (2)  is  a  first  order  differential  equation  relating 
p  to  p.  The  solution  of  (2)  is  given  by  a  functional  expressing  p  in  terms  of  the  history  of 
p.  The  incompressible  limit  arises  by  considering  functions  <f>n  and  fn  such  that,  for  any 
given  p,  the  solution  of  (2)  tends  to  a  constant  p  =  po  as  n  — > -  oc.  For  example,  we  may 
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take  4>n(p)  —  Cn(p  —  po)/p,  equal  to  a  positive  constant,  and  let  C„  tend  to  infinity.  If 
Tj(p)  is  independent  of  n,  then  the  incompressible  limit  yields  the  classical  Navier-Stokes 
equation  with  constant  viscosity  r\  =  rj[po)-  We  may,  however,  consider  a  limiting  case 
where  not  only  <f^n/dp  becomes  infinite  as  n  — ♦  oo,  but  drin/dp,  and  perhaps  fn,  also 
become  infinite.  We  must  then  in  general  regard  r)  as  a  functional  of  the  history  of  p,  even 
in  the  incompressible  case.  In  this  paper,  we  restrict  attention  to  the  case  where  tj  is  just  a 
function  of  the  instantaneous  value  of  p.  Of  course  we  shall  assume  throughout  the  paper 
that  the  function  p(p)  is  always  positive. 

For  the  classical  Navier-Stokes  equations,  it  is  well  known  that  the  stationary  problem 
is  elliptic  in  the  sense  of  Agmon,  Douglis  and  Nirenberg  [l],  that  the  time-dependent 
problem  is  “parabolic”  and  that  both  Dirichlet  and  traction  conditions  on  the  boundary 
satisfy  the  appropriate  complementing  conditions.  This  is  of  fundamental  importance 
for  studying  the  regularity  of  solutions  to  the  stationary  problem  as  well  as  for  studying 
the  time-dependent  initial-boundary  value  problem.  If  the  viscosity  is  dependent  on  the 
pressure,  ellipticity  can  fail.  This,  and  the  complementing  condition  for  Dirichlet  and 
traction  boundary  conditions  will  be  discussed  in  section  2. 

If  ellipticity  fails,  we  can  not  expect  well-posedness  of  the  initial-boundary  value  prob¬ 
lem.  In  section  3,  we  will,  however,  show  a  local  existence  result  for  the  Dirichlet  initial- 
boundary  value  problem  under  the  assumption  that  the  initial  data  are  such  that  ellipticity 
is  satisfied.  The  main  complication  compared  to  ordinary  Navier-Stokes  theory  is  that  we 
can  not  use  the  usual  Hodge  projection  to  eliminate  the  pressure.  Instead,  we  must  solve  a 
nonlinear  elliptic  equation  for  eliminating  the  pressure.  The  linearization  of  this  equation 
at  the  current  velocity  and  pressure  defines  a  projection  operator  that  depends  on  the  ve¬ 
locity  field  and  the  pressure  and  assumes  the  role  of  the  Hodge  projection.  This  makes  the 
theory  considerably  more  complicated  than  that  of  the  ordinary  Navier-Stokes  equations. 
Another  difference  is  that  we  can  of  course  no  longer  add  a  constant  to  the  pressure  with¬ 
out  altering  the  flow.  However,  this  indeterminacy  in  the  classical  Navier-Stokes  theory 
has  an  analogue  in  the  general  case,  and  in  order  to  obtain  a  unique  solution,  we  must  for 
example  prescribe  the  average  pressure  as  a  function  of  time,  in  addition  to  the  initial  and 
boundary  conditions  for  the  velocity  field. 


2.  Ellipticity  and  the  Complementing  Condition 
The  Navier-Stokes  equations  have  the  form 

p{u  +  ( u  ■  V)u)  =  div  [r?(p)(Vu  +  (Vu)T)j  -  Vp+  /,  (3) 

div  u  -  0. 


We  can  rewrite  this  in  the  form 


p(u  +  (u  ■  V)li)  -  T](p)Au  -r  {^'(p)  Vu  +  (Vu)T)  -  l}Vp  +  /, 


(4) 


The  ellipticity  and  complementing  conditions  concern  the  terms  of  highest  differential 
order  in  the  equation,  and  the  coefficients  of  these  terms  are  treated  as  though  they  were 
constant.  That  is,  in  order  to  discuss  these  conditions,  we  have  to  look  at  the  problem 


pti  =  fjAu  -f  (A  -  1  )Vp, 
div  u  —  0. 


(5) 


Here  the  constant  matrix  A  replaces  »j'(p)(\7u  +  (Vu)rj,  and  it  is  therefore  symmetric  and 
traceless.  For  stationary  problems  we  have  u  =  0.  The  stationary  Navier-Stokes  system  is 
elliptic  if  equation  (5)  (with  u  =  0)  has  no  non-constant  periodic  solutions  on  all  of  space. 

It  is  not  difficult  to  see  that  this  is  the  case  if  and  only  if  A  —  1  is  negative  definite  (for 
the  ordinary  Navier-Stokes  equations,  A  =  0,  and  the  equations  are  always  elliptic).  To 
see  this,  let  us  first  assume  that  A  -  1  is  negative  definite.  We  multiply  the  first  equation 
of  (5)  by  (1  -  A)- *u  and  integrate  over  one  period.  After  integrating  by  parts,  this  yields 


I  £(d,».)(l  -  A), dz  =  0,  (6) 


which  immediately  yields  Vu  =  0.  For  the  time-dependent  problem,  we  obtain  in  the  same 
fashion  that 

u,(l  -  A )~luj  dx  =  -  /  ^2  -  A )~k'{djUk)  dx,  (7) 

and  hence  periodic  solutions  with  non-zero  wave  numbers  decay  exponentially  with  a  rate 
proportional  to  the  square  of  the  wave  vector,  i.e.  the  equation  is  “parabolic”. 

If  A  -  1  is  not  definite,  the  behavior  changes.  By  taking  the  divergence  in  the  first 
equation  of  (5),  we  see  that  the  equation  for  p  is  given  by 


i's/S 


div  j(A  -  l)Vp]  =  0, 


(8) 


and  there  are  now  nonconstant  spatially  periodic  solutions  to  this  equation.  As  a  con¬ 
sequence,  we  can  find  nonconstant  periodic  solutions  to  (5)  in  the  stationary  case,  and 
exponentially  growing  solutions  in  the  time-dependent  case.  Solutions  to  the  initial  value 
problem  are  not  unique,  and  if  inhomogeneous  terms  are  added  to  the  equation,  they  may 
not  exist.  We  conclude  that,  in  contrast  to  the  ordinary  Navier-Stokes  equations,  the  initial 
value  problem  is  not  always  well-posed.  Only  as  long  as  the  eigenvalues  of  the  symmetric 
part  of  the  velocity  gradient  are  less  than  can  we  exPect  t°  Prove  a  local  existence 

result.  If  this  condition  is  violated,  problems  of  nonexistence  and  nonuniqueness  occur 
in  the  constant  coefficient  problem  and  hence  can  be  expected  in  the  full  Navier-Stokes 
system  as  well. 

Roughly  speaking,  ellipticity  is  a  condition  concerning  the  local  behavior  of  solutions 
near  points  away  from  the  boundary  of  the  domain  in  which  the  equations  are  studied. 
Near  the  boundary,  however,  the  nature  of  the  boundary  conditions  is  important,  and 
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“well-behaved”  problems  are  characterized  by  a  restriction  on  the  boundary  conditions, 
in  addition  to  eliipticity.  This  extra  condition  is  known  as  the  complementing  condition. 
The  boundary  is  regarded  as  locally  flat,  and  one  looks  at  a  half-space  problem.  Again  we 
consider  the  terms  of  highest  differential  order  in  the  equations,  and  we  look  for  solutions 
of  the  half-space  problem  that  are  nonconstant  and  periodic  in  the  directions  parallel  to 
the  boundary  of  the  half  space  and  decay  exponentially  away  from  the  boundary.  The 
complementing  condition  for  the  stationary  problem  means  that  no  nontrivial  solutions 
of  this  kind  exist,  for  the  time-dependent  problem  the  appropriate  condition  is  again  the 
exponential  decay  of  such  solutions. 

Two  types  of  boundary  conditions  for  the  Navier-Stokes  equations  are  particularly 
important.  The  most  commonly  used  are  Dirichlet  conditions,  i.e.  u  =  0  on  the  boundary. 
For  these  conditions  the  complementing  conditions  are  trivial,  since  the  energy  estimates 
(6)  and  (7)  carry  over  without  any  change  to  the  half-space  problem.  For  problems  with  free 
surfaces,  traction  boundary  conditions  are  important.  The  discussion  of  these  conditions 
is  much  more  complicated,  and  we  shall  only  consider  the  two-dimensional  case. 

On  the  half-plane  x  >  0  of  i?2,  we  consider  equation  (5)  with  the  linearized  traction 
boundary  conditions 

.duj  du2.  .  A 

’’('5iT+ aT)  +  /,'2,,  =  0' 

2i,^  +  (A,,  -  l)p  =  0. 

We  use  Fourier  transform  in  the  y-direction  and  Laplace  transform  in  time  and 
solutions  proportional  to  eertetay.  We  then  have  the  following  set  of  equations 

poui  =  -  a2«i)  +  (An  -  l)p'  +  Al2iap, 

pou2  =  ??(u2  -  a2U2)  +  Ai2p'  +  (A22  -  1)*«P,  (10) 

u  j  +  iau2  =  0. 

The  boundary  conditions  at  x  =  0  read 

r}(iaui  +  u'2)  +  Ai2p  =  0,  (11) 

2>7ui  +  (An  —  l)p  =  0. 

Here  '  denotes  differentiation  with  respect  to  x.  The  complementing  condition  says  that 
for  Re  a  >  0  and  0  /  0  there  are  no  nontrivial  solutions  which  decay  as  1  -*  00. 

Without  loss  of  generality,  we  can  assume  that  rj  =  p  =  a  =  1  (this  can  always  be 
achieved  by  rescaling).  Let  us  also  recall  that  An  +  A22  =  0,  and  that  eliipticity  implies 
that  Aj2  <  (1  —  An)(l  —  A22).  In  particular  this  means  that  An  and  A22  He  between  -1 

and  +1.  We  differentiate  the  first  equation  of  (10)  with  respect  to  x,  multiply  the  second 

by  ia,  and  add  them.  Using  the  divergence  condition,  we  obtain 


(9) 

look  for 


(An  _  1  )p"  +  2Ai2*p>  —  (A22  -  l)p  -  0. 


(12) 


The  solutions  to  this  equation  are  e~Xxt  where 


A  = 


-Aiji  ±  y/{An  -  l)(^22  -  1)  -  A22 

1  -  Au 


(13) 


Since  we  are  interested  in  solutions  that  decay  for  x  -*  oo,  we  must  choose  the  plus  sign 
in  the  numerator  of  (13). 

It  is  easy  to  see  that  p  =  0  does  not  lead  to  any  nontrivial  solutions  of  (10), (11),  hence 
we  set  p  =  e~Xx  with  X  given  by  (13).  From  (10)  we  find 


.  -JT+dx  ,  Un  ~  l)X-An%  _ Al 
Ul-Cie  +  -———e  , 


(14) 


_  „  - Wl+az  , 

uj  =  e2e  v  + 


(1  -  ^22)1  +  ->i 

A2 -l-o  6 


We  have  assumed  here  that  A2  ^  1  +  a.  The  case  A2  =  1  +  a  requires  a  separate  discussion, 
which  will  be  given  below.  From  the  incompressibility  and  boundary  conditions,  we  find 
three  constraints  for  the  two  parameters  Ci  and  C2.  This  leads  to  an  eigenvalue  problem 
for  0.  The  three  conditions  are 


— V^l  +  oc\  -f  tc2  —  0, 


ic  i  —  v/1  -t-  0C2  +  * 


(-4,11  ~  1)^  -  Ani 


A2  -  1 


A2  -  1  -  o 


(15) 


-2v^c,  -  2* -“.f  ~  +  (A„  -  1)  =  0. 

Let  now  7  =  %/l  4-  a.  After  some  algebra  equations  (15)  lead  to  the  following  equation  of 
third  degree  for  7 

(1  +  nr2) (A  +  ^(1  -  A„)  -  2* A  12(1  -  i2)  -  2A(1  -  Au)  +  2^(1  -  A22)  =  0.  (16) 

The  roots  of  (16)  were  computed  at  a  grid  of  points  in  the  set  {(An, -4 12) I  -  1  <  An  < 
1,A^2  <  (1  -  A fj)}.  No  roots  such  that  Re  -y  >  0  and  Re  a  =  Re  (i2  -  1)  >  0  were  found. 
If  we  assume  A  =  \/l  +  o,  we  arrive,  after  some  algebra,  at  the  equation 

i(A3  +  A)(l  -  An)  +  Ai 2(1  -  A2)  +  (An  -  Aa)iX  =  0.  (17) 


By  inserting  (13)  into  this,  and  comparing  the  real  and  imaginary  parts,  we  find  that  this 
is  the  case  only  if 

A 12  =  0,  and  A22  =  1.  (18) 


This  leads  to 


A  =  0,  0  -  -1, 


(19) 


and  ellipticity  does  not  hold  for  these  parameters.  We  conclude  that,  at  least  in  two 
dimensions,  the  complementing  condition  is  satisfied  as  long  as  ellipticity  holds. 


3.  The  initial  value  problem  for  Dirichlet  type  boundary  conditions 


We  consider  equation  (4)  on  a  bounded  domain  fl  C  Rs  with  smooth  boundary.  We 
want  to  show  the  existence  of  solutions  satisfying  a  given  initial  condition 


u(x,0)  =  t>(z), 


(20) 


and  Dirichlet  boundary  conditions 


u(x,t)  =  w(x,t)  on  dD. 


(21) 


It  is  assumed  that  (20)  and  (21)  are  compatible  with  the  incompressibility  condition  and 
with  each  other.  In  addition,  we  will  have  to  prescribe  the  average  pressure  as  a  function 
of  time, 

p(x,t)  dx  =  i{t).  (22) 


L 


n 

The  latter  condition  is  needed  to  remove  an  indeterminacy  already  present  in  the  ordinary 
Navier-Stokes  equations.  There  one  can  add  an  arbitrary  constant  to  the  pressure  without 
changing  the  flow.  If  we  want  to  make  the  pressure  uniquely  determined,  a  condition  such 
as  (22)  is  needed.  In  the  present  situation,  the  choice  of  t( t )  will  affect  the  velocity  as  well 
as  the  pressure. 

One  of  the  crucial  steps  in  dealing  with  the  Navier-Stokes  equations  is  the  elimination 
of  the  pressure.  In  the  classical  case  this  is  achieved  by  the  Hodge  projection  or,  in 
other  words,  by  solving  a  Neumann  boundary  value  problem  (see  e.g.  [9]).  In  the  more 
complicated  situation  of  equation  (4),  this  Neumann  problem  is  replaced  by  a  nonlinear 
elliptic  equation.  We  obtain  this  equation  by  taking  the  divergence  of  (4)  and  we  obtain  the 
associated  boundary  condition  by  multiplying  (4)  by  the  outer  unit  normal  of  the  domain. 
This  leads  to  the  problem 


div  jp(p)Au  +  (2p'(p)D  -  l)Vp}  =  div  jp(u  •  V)u  -  /}, 
with  the  boundary  condition 

n  •  jp(p)Au  +  (2rj'(p)D  -  l)Vpj  =  n  •  jp(u>  +  (u  •  V)u)  -  /}• 


(23) 


(24) 


We  want  to  solve  (23),  (24)  for  p  when  u  is  given.  We  require  that  p  satisfy  the  additional 
constraint 

[  p(x)  dx  =  -y,  (25) 

Ju 

where  7  is  a  given  real  number.  We  shall  prove  the  following  existence  and  uniqueness 
result. 


Theorem  1: 


Let  0  <  a  <  1.  Assume  that  u  is  in  H4(Cl)  D  C3,a(fl),  f  is  in  H2(Cl)  nC1,a(n),  w 
satisfies  the  condition  JgQw-n  dS  =  0  and  can  be  extended  to  a  divergence-free  vector  field, 
denoted  again  by  w,  which  lies  in  H2(d)  n  C1,a(0).  Moreover,  assume  that  p  is  a  smooth 
function  of  p  such  that  ri'(p)  is  bounded  and  Jimp—ioo  2^1  =  0.  Assume,  moreover,  that 
the  eigenvalues  of  D  are  strictly  less  than  1/(2  max r/),  if  maxr /  is  positive,  and  strictly 
greater  than  1/(2 min »/),  if  mini/  is  negative.  Then  (2S)-(25)  has  a  unique  solution 
p  £  H3(Q)  n  C2,a(Cl). 

Proof: 

We  set  g  =  p(u  •  V)u  +  pw  —  f.  To  find  solutions  of  (23)-(25),  we  look  at  the  family 
of  problems 

div  {rp(p)Au  +  (2ri/(p)D  -  l)Vp}  =  div  g,  (26) 

with  boundary  condition 

n  ■  {ttj(p)Au  +  (2rr/(p)D  -  l)Vp}  =  n  •  g,  (27) 

and  the  constraint  (25). 

For  r  =  0,  we  have  the  ordinary  Neumann  problem,  and  we  have  a  unique  solution 
p  6  H3(U)nC2'a(Cl)  for  every  given  g  £  H 2(fl)nC1,a(n)  and  76/?.  As  t  is  increased,  this 
solution  can  be  uniquely  continued  to  a  family  of  solutions  p,,  unless  one  of  the  following 
happens. 

1.  Ellipticity  is  lost. 

2.  The  norm  of  pr  in  H 3  n  C2,a  becomes  unbounded  as  r  approaches  some  value  rc. 

3.  The  linearized  problem  has  nontrivial  solutions. 

If  none  of  these  three  happens  up  to  r  =  1,  a  soluxion  of  (23)-(25)  exists.  Conversely, 
if  every  solution  for  r  =  1  can  be  continued  back  to  t  =  0,  we  have  a  uniqueness  result. 

Ellipticity  means  that  the  matrix  2ri/(p)D  -  1  remains  negative  definite.  According 
to  our  assumption  on  the  eigenvalues  of  D,  ellipticity  holds  up  to  t  =  1. 

The  estimates  of  Agmon,  Douglis  and  Nirenberg  can  be  used  to  show  that  the  norm 
of  p  in  H3  n  C2,a  remains  bounded  as  long  as  the  C'-norm  remains  bounded  and  the 
ellipticity  constant  remains  strictly  positive.  We  can  then  use  theorem  2.1  on  page  476  of 
1 4 j  to  obtain  a  bound  on  the  C^-norm  provided  there  is  a  bound  on  the  maximum  norm.  If, 
moreover,  we  assume  that  i/  is  bounded,  then  lemma  3.1  of  [5]  applies,  and  we  can  replace 
the  maximum  norm  by  the  //-norm.  We  can  derive  an  a  priori  bound  on  the  L1- norm 
(in  fact  on  a  stronger  norm)  of  p  as  follows.  We  multiply  (25)  by  p  and  integrate  over  the 
domain  fi.  After  an  integration  by  parts,  we  have 

/  Vp(l  -  2r»/(p)D)Vp  -  rAup(p)Vp  dx  —  /  p  div  g  dx  -  /  pg-ndS.  (28) 

hi  Ju  Jdu 

According  to  our  assumption  on  17,  we  can,  for  any  t  >  0,  find  a  constant  C(e)  such  that 
r?(p)  <  C(t)  -  <|p  .  When  we  use  this  estimate  in  the  second  term  on  the  left  of  (28),  we 
get  an  estimate  on  the  //’-norm  of  p. 


Finally,  we  have  to  show  that  the  linearized  problem  can  have  only  the  trivial  solution, 
i.e.  that  there  can  not  be  nontrivial  solutions  <f>  of 

div  {»V(pr)Au^  +  2Tri"(pr)DVpT<f>  +  (2rij'(pT)D  -  1)V<£}  =  0,  (29) 

with  boundary  condition 

«  •  {iV(Pr)Au<£  +  2TT]"(pT)DVpT<f>  +  (2rf?'(pr)D  -  l)V<j>}  =  0,  (30) 

and 

I  4>  dx  =  0.  (31) 

J  n 


This  problem  is  of  the  form 
with  boundary  condition 


div  {AV<£  +  f<t> }  =  0, 


(32) 


n  •  (AV<£  +  f<f>}  —  0,  (33) 

where  A  is  a  positive  definite  symmetric  matrix.  Let  us  consider  the  parabolic  equation 


4>t  =  div  (A V<£  +  f<j> }  (34) 

with  boundary  condition  (33)  and  let  T  be  the  time  1  evolution  operator,  i.e.  the  operator 
that  takes  an  initial  datum  <f>(t  =  0)  to  <f>(t  =  1).  It  is  easy  to  verify  that  Jn  T<f>  dx  = 
fa  <P  dx.  It  follows  from  the  strong  maximum  principle  for  parabolic  equations  (see  chapter 
9  of  [6])  that  T<t>  >  0if<£>0  and  <f>  does  not  vanish  identically.  This  implies  that,  if 
T<t>  —  <j>,  then  <f>  must  be  either  positive  or  negative.  Assume  the  contrary  and  let  <f>+  be 
the  positive  part  of  < <f>.  If  <f>  is  not  positive  or  negative,  then  <t>+  is  not  identical  to  <f>  and 
not  identical  to  0.  It  follows  that  T<f>+  >  T<f>  =  <f>  and  T<{>+  >  0,  hence  T<)>+  >  <f>+.  But 
this  is  a  contradiction  since  T<j>+  must  have  the  same  integral  as  <j>+.  Hence  the  equation 
T<t>  —  <f>  has  at  most  one  linearly  independent  solution.  Since  T  is  a  compact  operator  in 
L2,  the  operator  1  —  T  has  index  zero,  and  since  it  maps  to  functions  with  zero  integral, 
it  is  obviously  not  surjective.  Hence  the  nullspace  of  T  is  one-dimensional.  The  solutions 
to  (32),  (33)  satisfy  T<t>  ~  <j>,  and  can  therefore  not  have  zero  integral.  This  completes  the 
proof. 


Unfortuna‘ely,  the  assumptions  on  the  function  r?(p)  in  theorem  1  are  not  very  realistic 
from  a  physical  point  of  view.  The  experiments  [3]  show  that  rf'  increases  with  p.  Hence 
it  may  not  be  realistic  to  assume  that  r)'  is  bounded,  and  it  is  certainly  not  realistic  that 
v{p)/p  tends  to  zero  as  p  — *  +oo.  Negative  pressures  lead  to  cavitation  in  reality,  and 
hence  any  assumptions  on  r?(p)  for  p  — *  -oo  are  academic.  Under  realistic  assumptions, 
we  can  therefore  not  guarantee  the  existence  of  a  pressure  for  a  prescribed  velocity  field 
and  body  force.  Of  course  an  implicit  function  argument  can  be  used  if  the  velocity  and 
body  force  are  small. 

We  now  turn  to  the  discussion  of  the  Dirichlet  initial  value  problem  given  by  (3)  and 
(20)-(22).  We  make  the  following  assumptions  on  the  data  of  the  problem.  All  these 
assumptions  are  assumed  to  hold  on  some  time  interval  j0,<o]- 
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(i)  The  body  force  /  is  a  uniformly  Holder  continuous  function  of  time  taking  values 
in  /f2(ft).  Its  time  derivative  is  uniformly  Holder  continuous  with  values  in  L2(fl). 
Moreover,  at  t  =  0,  /  lies  in  C1,a(f2). 

(ii)  The  function  w  satisfies  the  condition  /an  w  •  n  dS  =  0.  It  can  hence  be  extended 
to  a  divergence-free  vector  field  in  the  interior  of  fl,  which  we  denote  again  by  w. 
We  assume  that  w  is  uniformly  Holder  continuous  with  values  in  H4( H),  its  time 
derivative  is  uniformly  Holder  continuous  with  values  in  if 2 (ft)  and  its  second  time 
derivative  is  uniformly  Holder  continuous  with  values  in  L2{U).  At  t  =  0,  w  lies  in 
C3'“(n).  Moreover,  w  lies  in  C1,OI(n)  for  t  =  0. 

(iii)  v  is  divergence  free  and  lies  in  H4{ f!)  n  C3,a  (n). 

(iv)  The  function  7  has  a  uniformly  Holder  continuous  derivative. 

(v)  For  t  =  0,  equations  (23)-(25)  (with  u  =  v)  have  a  unique  solution  po  €  H3{ fl)  n 
C2,a(n)  (e.g.  if  the  assumptions  of  theorem  1  hold).  The  matrix27?'(p0)[Vt;+(Vt;)T]- 
1  is  strictly  negative  definite  on  all  of  fl. 

(vi)  v(x)  =  tu(x,0)  on  an,  and  u(x,t  =  0)  as  computed  from  (3)  agrees  with  ti>(x,0)  on 

an. 

Theorem  2: 

Assume  that  assumptions  (i)-(vi)  hold.  Then  equation  (S)  with  the  data  (20)-(28) 
has  a  unique  solution  on  some  time  interval  [0,  Ti)  with  0  <  T\  <  to-  The  solution 

is  such  that  u  e  c°((o,r1),Jf4(n))  n  c1([o,r,),/f2(n))  n  c2([o,r,),L2(n))  and  P  e 
c°([o,r1),/f8(n))nc1((olr1),ffl(n)). 

Proof: 

We  shall  look  not  only  at  equation  (3)  itself,  but  also  at  the  time  differentiated  equa¬ 
tion.  Let  us  denote  a  =  u  -  Xu  (A  is  a  constant  to  be  specified  later)  and  q  =  p.  On 
differentiating  (3),  we  obtain  the  following  evolution  problem  for  a 

p\a  +  \a  4-  A2u  +  ((a  +  Xu)  ■  V)u  4-  (u  •  V)(a  4-  Au)]  = 

div  [»?'(p)?(Vu  +  (Vu)T)  4-  rf(p)(V(a  +  Au)  +  (V(a  +  Au))r)]  -Vq  +  },  (35) 

div  a  =  0. 


Mr  .. 


\\\.* 


Moreover,  we  have  the  boundary  condition 


a{x,t)  =  u>(x,t)  -  Au;(x,/)  on  an, 


and  the  average  of  q  is  given  by 


j  q(x,t)  dx  =  7 (t). 
J  n 


Finally,  an  initial  condition  for  a  is  obtained  from  equation  (3) 

a(x,0)  =  -Av  -  (v  •  V)v  +  l-  {div  ]*?(p0)(Vv  4  (Vr)r'  -  Vp0  -  /(x.0)}. 


V>\V. 

m 


it 


>■ . 


.v 

V-V-V-V -  -  A". 


Here  po  is  the  initial  pressure  which  is  guaranteed  to  exist  by  assumption  (v). 

The  strategy  is  now  to  express  u  and  p  in  terms  of  a  using  equations  (3),  (21)  and 
(22).  By  inserting  the  result  into  (35),  we  get  a  problem  involving  only  a  and  q.  From  this, 
we  eliminate  q  by  solving  a  Neumann  problem  as  above.  This  finally  yields  an  evolution 
equation  for  a,  which  we  shall  treat  using  the  theory  of  analytic  semigroups. 

At  time  t  =  0,  we  have  (3),  (21)  and  (22)  satisfied  with  u  =  v,  p  =  po  and  a  =  a(x,0) 
as  given  by  (38).  In  the  neighborhood  of  these  data,  we  want  to  show  that  if  A  £  R  is 
chosen  large  enough,  then  we  can  uniquely  resolve  the  equations  for  u  £  if  4  and  p  6  H3 
in  terms  of  a  €  H2,  /  £  if2,  w  £  H4  and  7.  In  a  weaker  norm,  we  shall  get  u  £  H3 
and  p  £  H 2  in  terms  of  a  £  if1,  /  £  if1,  tv  £  H 3,  and  7.  In  order  to  show  this,  we 
use  the  implicit  function  theorem.  For  this,  we  have  to  look  at  the  linearized  equations. 
Linearizing  (3)  with  respect  to  u  and  p,  we  arrive  at  the  following  problem 

div  [q(po)(^«  +  (Vu)T)  +  rj,(p0)p(Vu  +  (Vu)T)J-Vp-p((u-V)u+(u-V)u)-pAu  =  g,  (39) 

div  u  =  0. 

We  have  to  show  that  (39)  together  with  the  inhomogenous  Dirichlet  conditions 


u  —  xv  on  <9f2, 

(40) 

f  pdx  =  7, 

(41) 

n 


has  a  unique  solution  u  £  H 4,  p  £  if3  for  every  g  £  H2,  every  tv  that  is  the  trace  of  a 
divergence-free  vectorfield  in  H4,  and  every  7  £  R  (plus  the  corresponding  assertion  for 
the  weaker  norms).  The  inhomogeneous  boundary  terms  can  be  absorbed  into  g  by  letting 
u  =  tv  +  u,  where  w  is  an  extension  of  the  boundary  data,  so  we  can  henceforth  assume 
that  tv  -  0,  and  also  that  7  =  0. 

We  shall  derive  estimates  that  imply  uniqueness  of  solutions  to  (39)-(4l)  for  large 
enough  A.  Moreover,  these  estimates  are  such  that  they  hold  uniformly  for  a  family  of 
problems  which  continuously  connect  (39)  to  the  standard  Navier-Stokes  equations  (we 
can  use  a  homotopy  that  takes  rj  to  a  constant  function).  The  results  of  chapter  12  in  [l] 
then  imply  existence. 

Hence  we  now  let  g  =  0  in  (39),  w  =  0  in  (40),  7  =  0  in  (41)  and  we  assume  that  A  is 
positive  and  large.  We  first  multiply  (39)  by  (l  -  r?,(po)(Vu  +  (Vv)T)}_1u  and  integrate 
over  the  domain  fi.  After  an  integration  by  parts  this  yields  an  estimate  of  the  form 
(provided  A  is  large  enough) 


I;«!Ih'  •+  <Ci!i p\Lz- 

The  elliptic  estimates  of  Agmon,  Douglis  and  Nirenberg  jl]  lead  to 

IiwIIh2  +  IIpI'w  <  c^Aptti!^  +  lipiij,;)- 


(42) 


(43) 
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By  taking  the  two  together,  we  arrive  at 

+  IIpIIj/*  +  <  c3\\P\\L,.  (44) 

We  can  use  the  same  procedure  as  in  the  proof  of  theorem  1  to  eliminate  the  pressure,  and 
we  note  that  the  equation  for  the  pressure  is  independent  of  A.  This  gives  p  as  a  linear 
function  of  u  and  a  straightforward  energy  estimate  leads  to 

IIpIIh*  <  c*\\u\\h*  (45) 

holds.  Hence  the  mapping  u  €  H2  — *  p  G  L2  is  compact.  It  follows  that  for  every  c  >  0 
there  is  a  constant  C(e)  such  that 

IIpIU2  <  +  C(e)||ti||i(a  (46) 

(use  approximation  by  operators  of  finite  rank  and  the  fact  that  L2  is  dense  in  the  dual  of 
H2).  We  now  obtain  the  desired  uniqueness  result  by  inserting  (46)  into  (44)  and  choosing 
c  and  A  in  such  a  way  that  eC3  <  1  and  C(€)C3  <  A. 

Using  the  same  steps,  we  can  obtain  a  resolvent  estimate  for  complex  A.  Let  A  be 
in  any  sector  £  of  the  complex  plane  which  excludes  the  negative  half-axis  and  let  |A|  be 
sufficiently  large.  Moreover,  assume  that  tv  =  0,  7  =  0.  Then  the  solution  to  (39)-(4l) 
satisfies  an  estimate  of  the  form 

IMff*  +  llpllzr*  +  IAINU*  <  C||*|U»  (47) 

with  a  constant  C  that  does  not  depend  on  A. 

We  now  turn  to  equations  (35)-(37).  According  to  the  above  we  may  regard  u  and 
p  as  having  been  expressed  in  terms  of  a  and  the  data  of  the  problem.  Moreover,  we 
can  subtract  appropriate  reference  functions  from  a  and  q  such  that  (36)  and  (37)  are 
transformed  to  a  homogeneous  form.  Finally,  we  can  eliminate  q  by  using  the  procedure 
in  the  proof  of  theorem  1.  In  this  way,  we  obtain  an  evolution  problem  involving  only  a 
as  an  unknown.  To  this  problem  we  can  now  apply  the  theory  of  Sobolevskii  [7],  The 
essential  assumption  of  this  theory  is  that  the  terms  of  highest  differential  order  generate 
an  analytic  semigroup.  The  terms  of  highest  order  are  those  terms  on  the  right  hand  side  of 
(35)  which  involve  second  derivatives  of  a  and  first  derivatives  of  q.  These  are  of  precisely 
the  same  form  as  in  (39),  and  hence  the  needed  resolvent  estimate  follows  from  (47).  The 
remaining  assumptions  of  [7]  concerning  smoothness  of  the  nonlinearities  and  regularity  of 
the  data  are  easily  checked.  This  completes  the  proof. 
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